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Abstract. We have collected data for 69 Damped Lyman-a (DLA) 
systems, to investigate the chemical evolution of galaxies in the redshift 
interval 0.0 < z < 4.4. In doing that, we have adopted the most general 
approach used so far to correct for dust depletion. The best solution, 
obtained through minimization, gives as output parameters the global 
DLA metallicity and the dust-to-metals ratio. Clear evolution of the 
metallicity vs. redshift is found (99.99% significance level), with average 
values going from ~ 1/30 solar at z ^ 4.1 to ~ 3/5 solar at z ~ 0.5. 
We also find that the majority of DLAs (~ 60%) shows dust depletion 
patterns which most closely resemble that of the warm halo clouds in 
the Milky Way, and have dust-to-metals ratios very close to warm halo 
clouds. 



1. Introduction 

Many models of galaxy formation and evolution in recent years take as a refer- 
ence for the observational side the results coming from QSO absorption studies 
and in particular those objects that show a large HI column density, namely, 
the damped Lyman-a systems (DLAs) with Nhj ^ 10^'' atoms cm~^ (see for 
example Pei et al. 1999). DLAs have been widely studied both because they 
are believed to be the progenitors of present day galaxies and because the large 
HI column density allows one to probe the chemical enrichment problem. DLAs 
constitute so far the best laboratory where to directly measure the heavy ele- 
ment enrichment in a large interval of evolutionary phases, and to understand 
the processes of star formation and metal pollution of the Universe. However, 
this kind of investigation requires a careful consideration of the effects of dust 
depletion suffered by DLAs (Pei et al. 1991; Pettini et al. 1994). We present the 
analysis of a sample of 69 DLAs in the redshift range 0.0 < z < 4.4 to investigate 
their chemical state. We find that, after allowance for dust depletion corrections 
which are obtained with a very general approach, the DLA population clearly 
shows a metallicity-redshift evolution. 
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Table 1. Mean element abundances relative to hydrogen 

[X/H] [X/Fe] 



Element 


DLA 


Mn 


-1.26 ±0.51 


Zn 


-0.96 ± 0.44 


Cr 


-1.36 ±0.45 


Fe 


-1.61 ±0.48 


Ni 


-1.66 ±0.45 


Si 


-1.28 ±0.58 


S 


-1.37 ±0.42 



Warm Halo" SMC* DLA Warm Halo" SMC'' 

-0.60 ±0.13 -l.lOjio " -0.08 ±0.21 +0.04 +0-lllo;o6 

-0.64j:°:l4 +0.58 ±0.29 - +0.57^°;°® 

-0.51 + 0.13 -1.131";" +0.17 + 0.09 +0.13 +0.08t°;°4 

-0.64 + 0.06 -1.2ll',;;J^ _ _ _ 

-0.84 + 0.07 -l-73lo;ig -0.14 + 0.19 -0.21 -0.52^°;°^ 

-0.28 + 0.19 -0.57l°" +0.43 + 0.18 +0.36 +0.641°°^ 

-0.04 + 0.20 - ■ - - - ■ 



a) Savage & Sembach (1996); 6) Measured toward Sk 108, Welty et al. (1997). 



2. Chemical abundances and evolution 

We have collected data from the literature for a DLA sample, which includes 
69 objects. This sample represents the largest and most complete sample of 
DLAs for which measurements of HI and hcavy-clcmcnt column densities are 
available. The ions considered for abundance measurements are Fell, Znll, CrII, 
Sill, Mnll, SII, Nill. These ions are the dominant contributors to the abundances 
of the corresponding elements in HI clouds with high column densities, because 
they all have ionization potentials below 13.6 eV. 

In Table 1 we give the mean metal abundances relative to hydrogen and iron. 
They are presented with the customary definition [X/H]/)^^ = log(Xj/Yj)/)i^ — 
log(X/Y)Q, where Xj and Yj are the ion column densities of element X and Y. 
For comparison, the mean abundances for Warm Halo (WH) clouds (Savage and 
Sembach 1996) and the Small Magellanic Cloud (SMC, Welty et al. 1997) are 
also shown. We note that globally DLAs show [X/H] and [X/Fe] abundance 
ratios more similar to those of SMC and WH clouds, respectively. This suggests 
that metal abundances in DLAs are the result of chemical enrichment processes 
similar from the ones operating in the SMC and that the most common depletion 
pattern operating in DLAs is similar to the one observed in WH clouds. 

Indeed, to derive a complete picture of the DLA chemical state, one must 
correct for dust depletion effects. Since every element considered is affected 
by dust depletion differently, one must consider all measured species simultane- 
ously. In the Milky Way, a number of depletion patterns have been identified, 
showing highest depletions in dense disk clouds and lowest depletions in low 
density, warm halo clouds (Savage & Sembach 1996). We make a simplification 
assuming that the depletion patterns in DLAs may be reproduced by one of 
the four depletion patterns identified for the MW: Warm Halo, Warm Halo + 
Disk (WHD), Warm Disk (WD) and Cool Disk (CD) clouds (Savage & Sembach 
1996), thus modifying the dust-to-metals ratio to obtain the best match with 
the observations. By means of a minimization procedure we determine the 
best fit DLA metallicities and the dust-to-metals ratios. 

Fig. 1 shows the metallicity as a function of redshift. Filled symbols repre- 
sent DLAs with three or more measured elemental abundances for which it has 
been possible to obtain a proper best fit solution (37 DLAs). For the cases with 
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Figure 1. Total (gas+dust) metallicity compared to solar for 

69 DLAs as a function of redshift. Filled symbols represent the 37 
DLAs with 3 or more measured elements, open symbols are DLAs with 
two (solid bars, 16 objects) or one (dotted bars, 16 objects) element 
measured. The dashed line is the best fit in the case of linear corelation 
and that gives a slope dlog ZoLA/dz = —0.30 zt 0.06. 



only two elements observed, each fit has a zero formal error and, therefore, a 
reduced cannot be calculated; thus, the best fit is considered less significant 
(16 DLAs, empty symbols and solid error bars). Finally, for the cases where only 
one element is measured, we estimate the metallicity assuming a WH depletion 
pattern (16 DLAs, empty symbols and dotted error bars). The combination of 
the largest sample available (69 DLAs), a large redshift baseline (0.0 < z < 4.4) 
and a more accurate dust correction applied have led to the unambiguous detec- 
tion of the redshift evolution of metallicity in DLA galaxies, with mean values 
around 1/30 of solar at z ~ 4.1 to 3/5 of solar at z ~ 0.5. We found a significant 
linear correlation of metallicity with redshift (99.99% significance level) with a 
slope of —0.301 zb 0.056, which is clearly not consistent with a null gradient, 
indicating genuine evolution of DLA metallicity with redshift. In Fig. 1 we also 
show six boxes centered on the weighted averages over Az = 0.74 intervals and 
whose vertical widths mark the corresponding ±la weighted dispersion. In ad- 
dition, we note that the vertical dispersion of points greatly exceeds the total 
errors, indicating that, although all DLAs align along an average trend, there 
are real differences among individual objects in either initial conditions, or time 
of formation, or metal enrichment efficiency, or all of the above. 

Pei & Fall (1995) consider that a mean heavy-element abundance in the 
interstellar medium of DLA galaxies is given by the ratio of the total metal 
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content to the total gas content (i.e. Z oc Nfngfalsjf^j^Hij), which means 
that large HI DLAs dominate when calculating the global metallicity. This kind 
of analysis has been performed on a sample of DLAs using the Znll absorption 
line and a null result has been found for the evolution (Pettini et al. 1999), 
and it is not disproved if our sample is used. However, the lack of evident 
evolution in this case appears to be due to the fact that those DLAs with large HI 
column density are concentrated in the central redshift region (84% of DLAs with 
log Nhi > 20.7 are in the bin 1.4 < z < 3.5). In other words, we cannot exchide 
the redshift evolution of metallicity in high HI DLAs till when the considered 
sample will be fairly distributed in the considered redshift range. If there were 
a metallicity evolution also for high HI DLAs, dust obscuration could be too 
large at z < 1 to make DLAs detectable (a DLA with log Nhi = 21.0 and 
metallicity log Zjjla/Zq = —0.5 would extinct about 1.5 magnitude in the Lya). 
Moreover, this is where UV observations, required to detect Lya absorption, are 
very sensitive to any appreciable dust opacity. 

In 37 DLAs for which we have obtained a best fit solution, we have found 
that the majority ('~ 60%) show a WH-like depletion pattern, then ^ 30% and 
~ 10% have a WHD-like and WD-like depletion pattern, respectively. The 
predominance of WH-like conditions may be interpreted as evidence for a lower 
density in the absorbing clouds relative to clouds in the MW disk and/or to 
more efficient dust destruction processes operating in DLA galaxies, possibly 
due to more active star formation that may produce a more turbulent, dust 
hostile environment. In our sample there are 31 DLAs with measured Si and 
for which it has to perform minimization. For 8 over 31 systems, there is 
evidence of deviation of Si with respect the expectations, having a mean value 
of 0.28 lb 0.05. No correlation with redshift and metallicity is revealed. This 
result may be indicative of a-element enrichment in a fraction of DLAs, or of a 
lower efficiency of silicate formation than is found in the Milky Way. 

We like to stress that it is the combination of a large sample of DLAs 
(69 objects), a wide redshift baseline (0.0 < z < 4.4) and a more accurate 
dust correction applied to all heavy-elements measured, that has allowed us to 
demonstrate that the metallicity of DLAs evolves with redshift. The combined 
use of as many elements as possible allows one to much more efficiently track 
the metallicity evolution down to very metal poor gas clouds. 
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